Analytical computation of the magnetization probability density function for the 

harmonic 2D XY model 
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The probability density function (PDF) of some global average quantity plays a fundamental 
role in critical and highly correlated systems. We explicitly compute this quantity as a function of 
the magnetization for the two dimensional XY model in its harmonic approximation. Numerical 
simulations and perturbative results have shown a Gumbel-like shape of the PDF, in spite of the 
fact that the average magnetization is not an extreme variable. Our analytical result allows to test 
both perturbative analytical expansions and also numerical computations performed previously. 
Perfect agreement is found for the first moments of the PDF. Also for large volume and in the 
high temperature limit the distribution becomes Gaussian, as it should be. In the low temperature 
regime its numerical evaluation is compatible with a Gumbel distribution. 
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The probability density function (PDF) of some global 
average quantity plays a fundamental role in critical and 
highly correlated systems. For example, it has been used 
in the so-called hypothesis of universal fluctuations [l|, 
which states that the PDF for critical systems -like the 
two dimensional XY model- as a function of the cen- 
tered parameter of order divided by its standard devi- 
ation should be a universal function F of a Gumbcl- 
type. This claim should hold independent of the size and 
temperature for the magnetic systems and Reynold num- 
ber for turbulent experimental systems. Approximated 
Gumbel distributions appear in quite dissimilar physical 
systems, such as a confined turbulent flow or the rough- 
ness 1// noise in a resistor [2j]. 

There was some evidence against this hypothesis re- 
ported in In particular, in Ref. [i| for instance, 
an analytical expression for the PDF for the full 2D XY 
model was computed systematically by means of the loop 
expansion, and the validity of the " generalized universal- 
ity" has been linked to renormalization group (RG) prop- 
erties. The 2-loops analytical expression for the PDF 
shows an explicit temperature dependence. As a con- 
sequence, its skewness and kurtosis computed perturba- 
tively up to two loops (first order in T) show an explicit 
temperature dependence. More recently ||, the third 
and fourth normalized moments of the PDF -the skew- 
ness and kurtosis respectively- were computed analyti- 
cally in the two-dimensional harmonic XY model. Their 
explicit temperature dependence was explicitly demon- 
strated, which holds even in the thermodynamic limit. 
This result was an indirect analytical proof of the failure 
of the claim of universal fluctuations, allowing therefore 
to confirm the explicit temperature dependence of the 
PDF reported in 

In spite of this, further papers supporting the proposed 
" generalized universality" were published 6] , raising the 
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question whether there is a underlying mechanism re- 
sponsible for this approximated phenomenon. 

In a recent paper Q the shape of the different dis- 
tributions appearing in several critical systems has been 
phenomenologically linked to scaling arguments and to 
the concept from Renormalization Group Theory of clas- 
sification of scaling variables as irrelevant, marginal and 
relevant. Nevertheless, one should note at this point that 
the so-called "generalized universality" is not related to 
the concept of exact universality of RG, because it only 
holds approximately, as it was explicitly shown in ref. 
[H by using the loop -or temperature- expansion in the 
context of the two-dimensional XY-model. 

It is further stated in 0] that an explicit analytical 
computation of the PDF for the two dimensional XY 
model is still missing. This claim also motivates our 
present computation. 

Now we will explain how to compute a general ana- 
lytical expression for the PDF of the magnetization by 
using the spin wave or harmonic approximation of the 
2D XY model. This expression is valid for arbitrary 
system size L and temperature T. The 2D XY model 
consists of planar spins cj> x defined on a periodic two di- 
mensional square lattice A oi N — L 2 lattice sites, which 
are coupled with nearest neighbors by cosine interactions. 
According to RG arguments Q, in the low temperature 
phase and sufficiently below the Berezinskii-Kosterlitz- 
Thoulcss critical temperature, the physics of this model 
is entirely described by its harmonic approximation, the 
2D HXY model. Indeed, in [J] it has been explicitly 
shown -by means of the loop expansion- that the effect 
of the anharmonic corrections to the spin wave approx- 
imation on the PDF is merely a renormalization of the 
temperature. Nevertheless, and in spite of this fact, one 
should take into account the periodicity of the variables 
of the XY model for the boundary conditions, which 
leads to the contributions coming from winding configu- 
rations [9[ . This contribution turns out to be numerically 
very small and therefore one expects to obtain with the 
present model a trustable numerical approximation of the 



2 



2D XY model in the large volume limit and in the low 
temperature phase. 

The Hamiltonian of the 2D HXY-model is up to a con- 
stant 



H{cj>) = i J {<f>, -A<f>) 



(1) 



where A is the Laplace operator on the lattice, J is the 
ferromagnetic constant and (</>, ip) = ^2 0( x )*'/ 3 ( x ) stands 
for the scalar product on the lattice. We use a system of 
units where Boltzmann's constant is set equal to unity 
throughout the computations and identify T with the re- 
duced temperature T/J. Although this model has no 
phase transition, it is a critical model in the sense that it 
has an infinite correlation length. This Gaussian model 
has been analytical extensively studied, because it repre- 
sents the starting point for perturbative expansions and 
due to the involvement of Gaussian integrals. Some use- 
ful physical quantities can be expressed in terms of the 
Fourier representation of the lattice propagator G 



G(x) 



1 

TV 



E 

(Kl )V0 



exp(— iK • x) 
(Kl) 2 



(2) 



where is the lattice momentum defined as usual as 
{K^ji = 2sin(Ki/2), with i = 1,2 and Ki lies in the 
first Brillouin zone, Ki = (2n/L)n with n S 1 and 
— 7r < Ki < it. The sum runs over all possible values 
of Ki for which (K^) 2 does not vanish. This comes from 
the fact that the Goldstone mode, which is originated by 
the invariance of the original Hamiltonian under a global 
rotation of the spin variables, and which leads to the 
translation invariance on the lattice, must be removed 
from the calculation. As it was first shown in Ref. 
the PDF can be defined as the Fourier transform of the 
partition function Z(q) of an auxiliary theory, which dif- 
fers from the original theory by a dimension perturba- 
tion with a very small imaginary coefficients iq/N. This 
theory turns out to be asymptotically free in the infrared 
(this is the case of the 2D XY model), viz. 



P(M) 



dq 
2^ 



exp{iq(M - (M)}Z(q) (3) 



where 



Z(q) = 



,-iq(M) 
Z 



D(f>exp 



(4) 

with partition function Zq = J D4>exp{—[3H (</>)}. 
The mean (M) and the higher order moments of the 
PDF are obtained as usual as the integrals (M p ) = 
J M p P(M)dM . In particular, the mean square fluctu- 
ation is defined by a 2 — ({M — (M)) 2 ). The partition 
function Z(q) can be written as an integral over the nor- 
malized Gaussian measure d^TG(4>) of covariance TG: 



Z(q)=exp{-iq(M)) 



d[i TG {(j)) exp <{ — ^2 cos ' 



(5) 



where G denotes the lattice propagator given by eqn.([2]). 

In terms of the PDF, the hypothesis of BHP univer- 
sality stipulates that the PDF considered as a function 
of the reduced variable /i = (M — (M))/a is a univer- 
sal function F, and should be the same for a wide class 
of strongly correlated critical systems, independent of the 
system temperature T and volume N, provided N is large 
enough: 



F(jj) - oP 



(M - (M)) 



(G) 




FIG. 1. In this figure we plot the approximate expression 
for Z(q) given by equation (|15|l for L = 16. Z(q) is also 
evaluated numerically from eqn. (|10[l for a lattice with L = 16 
and T = 10. In order to compare both results and to obtain 
an analytical simple expression a Gaussian fit is performed 
obtaining x 2 — 1. 

Moreover, this universal function should be non- 
Gaussian and has a functional dependence very close to 
the one corresponding to a Gumbel distribution for cer- 
tain exponent a, which is known in the context of ex- 
tremum statistics [10j . We notice that up to the first fac- 
tor, the partition function Z(q) defined by the integral ^ 
corresponds to the partition function of the Sine-Gordon 
theory with a small imaginary fugacity z — iq/N . It is 
straightforward to see that Z(q) is the generating func- 
tional of the centered high order moments of P, i.e. 



d n 



Z(q) 



((M-(M)r) 



(7) 



g=0 



Motivated by the method used to study the renormal- 
ization group flow in the Sine-Gordon model fill ] . we 
will use the Fourier representation of an exponential fac- 
tor (which is also known from the high temperature -or 
cluster- expansion of the XY model [12|). 



exp{iA cos (j>} 



oo 

E 



z m J m (A)exp{im0}, (8) 
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where J m (A) are the Bessel function of integer order m, 
and the functional identity for thegenerating function 
of the Gaussian measure dnTG{4>) [13j | : 

d^Toisf) ex P (« (</>, /)) = exp(-i (/, TGf)), (9) 

to compute analytically the integral of eqn. J5J. The re- 
sulting expression for Z(q), which is our main analytical 
result, reads 



Z(q) = e 



- p -iq(M) 



N 



mi,m 2 ,..,mjv k=l 



exp < — T2_, mimjGfei — Xj) 



(10) 



From this explicit expression for Z(q) and using eqn. (JT]) 
one can compute explicitly the moments of the PDF, 
obtaining for instance 



(M) = exp 



TG(Q) \ 
2 J 



(M) 



zeA 



1/2 



(M 3 ) 



2JV 2 



E 

x,yGA 



(11) 



(12) 



(13) 



(e- TG W cosh{T[G(y)+G(x-y)]} 
+e TG Wcosh{T[G(y)-G(x- y )]} 



These expressions are exact and agree with previous re- 
sults reported in Q. Higher centered moments of the 
PDF can be computed as well by using equations ([7]) 
and (JTUf . The result perfectly agrees with their corre- 
sponding expressions reported for instance in Ref. @ 
(see eqns. (2.4) and (2.5)). 

In order to obtain the analytical expression for the 
" universal distribution function F" one has to insert the 
expression for Z(q) given by eqn. (|10p into ([3]), which 
leads to the result 



dQ 
2tt 



e^{iQfi}Z(-Q/a), (14) 



where (M) and a are given by eqns. (fTTj) and (p~2|) . 

Now we want to study numerically the expression for 
the PDF deduced for the 2D HXY model as a function 
of the system temperature T and volume N. From the 
equation (fTTj) it follows that in the high temperature 
limit, the mean magnetization goes to zero. In this limit, 
the exponential factor of equation pop vanishes for all 



values of TOfe € Z except for nik 
limit equation (flU)) goes into 



0. Therefore, in this 



Z(q)\ 



N 



(15) 



Moreover, for large volume the above expression for 
Z(q) becomes Gaussian. This point can be shown an- 
alytically performing a Taylor expansion valid for small 
arguments of Jo(q/N). Within this approximation, the 
mean square fluctuation of Z(q) is given by a w \^2N. 
On the other side, Z(q) can be also evaluated numeri- 
cally by performing the sums appearing in eqn. (|10p over 
the relevant configurations of m^-values. Both curves 
are displayed in Fig[TJ for a square lattice of lattice size 
L = 16 and T = 10. They felt onto the same curve with 
remarkable accuracy. Perfect agreement is found when 
a Gaussian distribution is fitted to Z(q), with a value of 
X 2 = 1. Also the numerical value for a — 22.51 agrees 
with the analytical expression for a in the high-T limit. 



F[Z(qfo)] for T=1 




FIG. 2. This figure shows the inverse Fourier transform of 
Z(q) -or equivalently the PDF defined by eq. (|14[) - corre- 
sponding to the data appearing in figure [T] 

The PDF itself can be computed both numerically and 
analytically in the high temperature limit. Indeed, using 
the approximated numerical values obtained already for 
Z(q), we evaluate its Fourier transformation defined by 
the integral of Eqn. ([3]). One may use instead the ac- 
curate analytical Gaussian expression for Z(q) found in 
this limit to perform analytically the integral of Eqn. ([3]). 
Both results are in perfect agreement as it is shown in 
figure El where a lattice of lattice size L = 16 and tem- 
perature T = 10 was used. The values obtained for the 
skewness and kurtosis, s « and c « 3 respectively, per- 
fectly agree, within only a few percent of error with their 
corresponding values reported in Ref. || . 

Finally, we can obtain the PDF itself in the low tem- 
perature regime using direct numerical integration. In 
figure 3 a numerical evaluation of the distribution for the 
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FIG. 3. This figure shows with a dashed line the numerical 
inverse fourier transform of Z(q) given by eq. (|14|l for lattice 
size L = 16 and temperature T = 0.7. For comparison we 
plot (full line) the analytical expression reported in ref. 
found for the roughness 1/f noise. The agreement between 
both curves is noteworthy. 



reduced magnetization is shown for a lattice of lattice size 
L = 16 and temperature T — 0.7 (dashed line). This cor- 
responds to the inverse fourier transform of Z(q) given by 
cq. (|14p . In order to compare our expression with the an- 
alytical expression for the roughness of signals displaying 
1/f power spectra reported in ref. [2| , we also plot (full 
line) the corresponding Gumbel-like distribution. The 
agreement between both curves is noteworthy. 
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